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ABSTRACT 


The axial corner flow is analyzed for the imcompressible laminar boundary 
layer flow. The governing equations are derived from the Navier-Stokes equa- 
tions by neglecting second derivative terms of the axial direction. An 
alternating direction implicit method is used to solve the equations in 
/ primitive variables. 

For large values of the axial distance, the present calculation for the 
case of coplanar leading edges is in good agreement with those of previous 
investigators . 
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INVESTIGATION OF AN INCOMPRESSIBLE FLOW ALONG A CORNER 
BY AN ALTERNATING DIRECTION IMPLICIT METHOD 


By 

Dashrath K. Patel ^ 


SUMMARY 

The axial corner flow is analyzed for the incompressible laminar boundary 
layer flow. The governing equations are derived from the Navier-Stokes equa- 
tions by neglecting second derivative terms of the axial direction. The- complete 
set of governing equations is prepared in primitive variables. In this first 
formulation the axial pressure gradient is assumed to be zero. An alternating 
direction implicit method is used to solve the governing equations. For rapid 
convergence simultaneous solution is obtained by tridiagonal block matrices. 

Results are presented for the case of coplanar leading edges in the first 
phase of computation. For large values of axial distance these results are in 
good agreement with those of previous investigators. For small values of axial 
distance, the present calculations indicate that the similarity concept used by 
previous investigators is not valid. 

INTRODUCTION 

The three-dimensional viscous interaction occurring in axial corners, 
formed along the interaction of two plane surfaces, is of fundamental interest 

in aerodynamic configurations such as wing-body and wing-fin jxmctions, wind 
tunnel .and its support, and other applications. 

The mutual interaction of the boundary layers on the two intersecting sur- 
faces must be taken into consideration, even in the lowest order approximation 
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for such a flovf. This makes it essential to consider the comer flow problem as 
a three-dimensional problem. 

The early analysis of Carrier (ref. 1) for the comer flow problem has been 
recently improved upon by Rubin (ref. 2) , who formulated it as a singular pertur- 
bation problem, including the three-dimensionality of the configuration. The 
method of matched asymptotic expansions was used to match systematically the 
corner' flow to the boundary layer for large lateral distances from the corner 
and to the potential flow for large distances away from the walls. Pal and 
Rubin (ref. 3) determined the asymptotic features of the corner flow and demon- 
strated that consistent asymptotic series exist for the "corner layer." Further, 
all flow variables were found to exhibit algebraic decay into the boundary 
layer away from the corner; also, the crossflow velocities were shown to decay 
algebraically into the potential flow. Numerical solutions for the incompressible 
comer flow were obtained by Rubin and Grossman (ref. 4) using the Gauss Siedel 
method of successive iterations and by Ghia and Davis (ref. 5) using the ADI 
method in similarity variables. 

These flows have certain characteristics which are common to boundary layer 
flows; for example, the velocities in the plane normal to the mainstream direction 
are usually much smaller than the main flow velocity. Also, one expects gradi- 
ents that exist in the cross directions to be larger than the gradients in the 
main flow directions. Therefore, it might be reasonable to use a boundary layer 
scaling on the Navier-Stokes equations to effect some simplification; however, 
this procedure yields an inconsistent set of equations. To lowest order, the 
cross -stream momentum equations reduce to statements that the pressure is 
constant, and the resulting niimber of equations is not sufficient to determine 
the remaining unknown quantities. A set of equations which is a higher order 
approximation to the Navier-Stokes than the boundary layer equations, the 
parabolic-elliptic Navier-Stokes (PENS) equations, is necessary to overcome 
these difficulties. 

The PENS three-dimensional equations are given in reference 6, where the 
leading-edge flow on a finite-width flat plate is calculated. The alternating 
direction implicit (ADI) method, introduced in reference 7, is used to solve 
the equations. The two- and three-dimensional problems are reviewed in reference 
8, where a new predictor-corrector method is presented as the most appropriate 
one to use. 
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The PENS concepts have been utilized in subsonic, incompressible flow for 
internal (ref. 9) and external (ref. 10) cases. However, the need to solve a 
Poisson equation for the cross-plane pressure in these incompressible flows 
separates them from the procedures used in references 6 and 8. 

The character of the equations is such that they are well suited to solution 
by the ADI method. This finite difference method generated block -tridiagonal 
matrix equations which allow for the simultaneous solution of all the variables 
at each location. 

This solution procedure of the PENS system is supplied to the corner layer 
problem. Computations are made in the cross-plane and the cross-planes are 
coupled by a marching procedure along the axial direction. 

SYMBOLS 

a^, . coefficients of matrix (A) in typical difference equation (37) 

^ defined in equation (40) 

b. . coefficients in matrix (B) defined in equation (37) 

J 

c. . coefficients of matrix (C) defined in equation (37) 

^ > J 

d. coefficient of vector D defined in equation (37) 

f(ri) Blasius function 

g(n) function defined in equation (26) 

i axial coordinate along x for discretized problem 

j lateral coordinate along n for discretized problem 

k lateral coordinate along ? for discretized problem 

L reference length 

M Mach number 

N value of j or k on the free stream end 

P static pressure 

Q dummy variable used for coefficients of derivatives 

Re Reynolds number 

u streamwise velocity in x direction 
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transverse velocity in y direction 
lateral velocity in z direction 

vector of four variables (u, v, w, p) , see equation (37) 

Cartesian coordinates 

step size along x direction, see equation (31) 
stretched transverse coordinate, defined in equation (7) 
stretched lateral coordinate, defined in equation (7) 
constant defined in equation (26) 
total boundary layer thickness 

parametel* in series expansion, defined in equation (6) 

similarity coordinate in lateral direction, defined in equation (16) 

step size along ? direction, defined in equation (31) 

similarity coordinates in transverse direction, defined in equation (16) 

step size along n direction, defined in equation (31) 

ratio of the consecutive step sizes defined in equation (31) 

shear coefficient of viscosity 

density 

constant- defined in equation (26) 

diommy variable used for u, v, w, and p 


Subscripts 


Q 1 2 variables for zero, first, and second order solutions of the flow, 

’ ’ ' see equation (6) 

“ free stream condition 

j denotes a grid point along n direction 

k denotes a grid point along ? direction 

max maximum value 

w wall value 
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step size behind grid point, i, j, k, see equation 
+ step size ahead of grid point, i, j, k, see equation (31) 

Superscripts 

i denotes a grid point along x direction 

n values of variable at iteration n, defined in equation (E-4) 

* . values of the variable at the intermediate step of ADI method. 

Also dimensionless corner flow variables defined in series 
solution in equation (6) . 

^ dimensional flow variables in (x,y,z) coordinates 

dimensionless flow variables in (x,y,z) coordinates 
-»• vector quantities 

' derivative of n or 5 

ANALYSIS 

The configuration of the viscous flow along a right angle corner formed by 
the intersection of two equal semi-infinite flat plates is shown schematically 
in figure 1. Following Rubin’s analysis (ref. 2), four different flow regimes 
are depicted. They are presented in figure 2. Region IV is primarily the region 
of present interest. 

Formulation of the problem for region IV requires knowledge of the following 
flows in regions I, II, and III. 

1 . The zero-order potential f loiv in region I 

2. The first-order boundary layer flow in regions II and III 

3. The first-order potential flow in region I, resulting from the 

normal flow at the outer edges of regions II and III. 

4. The second-order solution in regions II and’ III which is considered 
for crossflow velocities. 

A detailed asymptotic analysis of the corner layer equations is also needed 
to prescribe the proper far field boundary conditions for the "corner layer" 
flow. 
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Region IV depicted in figure 2 is referred to as the "comer layer" flow 
where the flow is truly three-dimensional ^ because of the interaction of the 
boundary layers in regions II and III. The governing differential equations 
for this comer layer flow will now be discussed. 

Governing Differential Equations 

The general dimensionless Navier-Stokes equations for a viscous, incompressible 
steady state fluid are 

Continuity 
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The dimensionless variables used in the above equations are defined as 
follows: 
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Now the corner layer flow solution is considered as the first terra in a 
series expansion, for the solution of the complete Navier-Stokes equations, 
using some appropriate small parameters e^. Therefore, the variables in 
equations ( 1 ) to ( 4 ) are represented by the following series; 

U = U* + El u| + £2 u|^ + ... 

^ = V* + El V? + S2 vJ + ... 

0 1 ^ 2 . 

■ (63 

W = W* + Si W* + £2 w| + ... 

P = Pq + Pf + £2 ?2 • • • 

where £i, £2 ... are the parameters of the series and £2 « £1 << 1 and 

so on (see ref. 11) . 

As discussed in the introduction, the cross velocities are usually much 
smaller than the main flow velocity. Therefore, a boundary layfer scaling is 
taken into account, and v* amd w* are assumed zero in equation (6) . 

With the boundary layer concept the independent lateral variables for the 
corner layer are stretched without change of axial length x follows: 

Y = y Re^/^ and Z = z Re^/^ (7) 


9 



Further, in equation (6) v* and W* are zero and thus v and W are of 
order in equation (6); From the continuit/ equation it can be shown that 


ei = Re- 1/2 

With consideration of the power series, it is assumed that 

^2 ~ ~ (9) 

This assumption is appropriate to obtain the governing equations by retaining 
the lowest order terms of the Navier-Stokes equations [eqs. (1) to (4)]. Thus, 
the following equations result: 

Continuity 


9uq 9vf 9wj 

"9F ~W Iz" 


0 


CIO) 


x-momentum 


u 


9uo 

0 


f\ * 

3uq 


* 

9uo 


+ ^ 

1 9Y 1 9Z 


9P0 

dX 




a^uo 


9Y2 3Z2 


Y- momentum 


3pS 


= 0 


3Pl 

'W 


= 0 


9Vi 9Vj 9Vi 

+ Y* ^ tj* 

0 9x 1 9Y 1 9Z 


9P2 


aY 


9Y^ 


az- 


( 11 ) 


( 12 ) 


(13) 
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Z -momentum 


3P0 

3Pl 


3Z ■ ^ 

3Z ■» 

(14) 

9wi 

3w^ 3W]^ 3p2 3^W| 


U* -T — + V* 
0 3x 1 

3Y ^ 3Z 3Z 3,2 

(15) 


Near the leading edge, the normal velocity is very large. Thus, the solution 
near the leading edge is singular in these variables. Therefore, to get the 
finite value of v and w near the leading edges, a transformation is applied 
to equations (10) to (15) . To get more details of the flow near the leading 
edge region, using a uniform or a variable grid in the computational plane, the 
lateral independent variables are also transformed. The system in (x,Y,Z) coor- 
dinates is transformed to (x,ri)?) coordinates as follows: 


u(x,n,^) = u*(x,Y,Z) 
v(x,n,?) = v|(x,Y,Z) 

(16) 

w(x,n,?) = m w*(x,Y,Z) 
p(x,n,?) = p|(x.Y,Z) 

where 

n = — S = — 


The governing differential equations for the comer layer in terms of the 
new coordinates (x,n,?), and the new dependent variables are written: 



Continuity 




3n 



9u 3w \ 
3? ■ 


0 


( 17 ) 


X-moraentum 


2 X U ^ - (UTl - V) ^ - (u? - w) 


3? 


l!iL 

3n^ 


3^u 

3?^ 


( 18 ) 


r|-momentum 


_ 3v ^ ^ 3v ^ ^ ^ 3v 

2X U ^ - (un - V) — - (us - w) — - uv 


^ 3p . 3^v . s^v 
= - 2 x -r^ + + 


3 ri 


3ti" 


9?' 


( 19 ) 


? -momentum 


o 3 w , . 3 w 

■ to - to 


f rr ^ 3W 

(u? - w) 9 ^ - UW 


- 3p 3^w 3^ 

= -2xt^ + + — 


9? 


9p^ 


■w 


9?^ 


( 20 ) 


There are four unknown variables, u, v, w, and p, and there are four 
governing equations. Thus, a solution is possible. But, pressure p appears 
only in first derivative and there is no p term in the continuity equation. 
Therefore, a numerical computation will create instability. To resolve this 
problem, a separate Poisson type pressure equation is derived as follows. The 
x-momentum equation is differentiated along x> tli® n-momentum equation is 
differentiated along n> the ^-momentum equation is differentiated along 
and then all three equations are added. With .the help of the continuity equation, 
the resulting equation is simplified, and finally the pressure equation is 
obtained as follows: 
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+ i!e. = 1 


3n^ 


8?^ 



8u 

3n 



3v 3w \ 
3n ^ 3? J 


3v 3w 3y 3w 
^ 3n 35 ~ 35 3n 


iH. 

3n 
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3v 3v 
3n “ ^ 35 
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Thus, the final form of the differential equations governing the corner layer 
flow is given by equations (18) through (21) . 

Boiindary Conditions 

As the flow equations (18) to (20) are of the PENS type and equation (21) is 
of the elliptic type in the cross plane (ri,5), the boundary conditions are required 
on all boundaries of a closed region in the (ti,5) plane. The boundary conditions 
along the plate surfaces n = 0, 5=0 are taken similar to no-slip boundary 
layer flow. 

On the far field boundary from the corner (in the boundary layer region) , 
the asymptotic analysis of the comer layer equations provides the necessary 
boundary conditions. For the solution of corner flow within a finite computational 
domain, the boundary conditions as provided by Pal and Rubin (ref. 3) or Ghia 
and Davis (ref. 5) will be used. The complete boundary conditions used for the 
domain are given below. 

At the plate surface, n = 0 (along OC in figure 2) 
u(0,5) = 0, v(0,5) = 0, w(0,5) = 0 ( 22 ) 


and 


3n 


(0,5) 


= 0 
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At the plate surface, ? = 0 (along OA in figure 2) 
u(n,0) = 0, v(ri,03 = 0, w(n,03 = 0 


( 23 ) 


and 




= 0 


At the far field boundary, the boundary conditions in boundary layer regitrn 
and in inviscid region are separately provided. It v-fas assumed that boundary 
layer edge occurs when — reaches a value of 0.99. 

. oo 

At the far field boundary, within the boundary layer region at ? = 
and for all n (along CD in figure 2) : 

u(ri,?) = f'(n) + X n f"(n) 


v(ti,c) = [n f>(n) - f(n)] 

- x[-o^ f”(n) + 3n f * (n) f(ri)]?”^ 

w(r),c) = e g(n) - 4x £'(n)r^ 


(24) 


and at n = n „ and for all ? (along AB ■ in figure 2) 

u(n,s) =£'(?) + X ? f"(c)n“^ 

w(n,c) = [? f’(?) - f(?)] 

- x[-?^ f"(S) + 3? f«(c) + f(?)]n"2 

v(n,e) = 3 g(0 - 4x f'(s)n“^ 

where 


(25) 
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£(n) is the Blasius function 
X = -2.5 


B = [n - f(riD] = 1.21678 ‘(26) 

g(Tl) = f'Cri) ? dx 

To calculate the boundary conditions for the inviscid corner region, the 
first-order potential flow solution in the region IV is essential. The solution 
for three velocity components is shown as follows: 
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The derivation and details o£ the above equations are represented in 
Appendix A. The pressure is calculated by the use o£ the Bernoulli's equation. 
The expression £or pressure calculation is given as follows: 



v^ 4- \ 

2Re X / 


C30) 


In the boundary layer region (along AB and CD in figure 2) , the presure 
is assumed constant. 

The mathematical model for the corner layer flow is now complete. The solu- 
tion of these nonlinear coupled partial differential equations [eqs. (18) to (21) 
subject to the boundary conditions given in equations (22) to (30), is obtained 
numerically using the ADI method. The details of the numerical analysis are 
presented in the next section. 

NUMERICAL ANALYSIS 
ADI Method 

Although explicit difference schemes are conceptually less complex and thus 
more easily coded than implicit methods, explicit finite difference methods are 
restricted to small axial distance steps relative to the spatial grid size for 
numerical stability. Consequently, explicit methods generally require a large 
total computation time to reach a stea.dy state flow condition. On the other 
hand, most implicit methods are found to be unconditionally stable in linearized 
stability analyses for model equations similar to the Navier-Stokes equations. 
Therefore, it is expected that implicit methods will allow larger axial steps 
than explicit methods when applied to the Navier-Stokes equations. 

The ADI technique developed by Peaceman and Rachford (ref. 7) was used in 
the present study. This ADI method is a two-step procedure requiring reduction 
of tridiagonal matrices for which an efficient solution algorithm, the Thomas 
algorithm, exists. Roach (ref. 12) states that ADI methods are currently the 
most popular approach to computing incompressible viscous Navier-Stokes flows. 

He notes that although these methods have not been found to be unconditionally 
stable in practice, the somewhat larger axial steps obtainable have resulted in 
faster overall computation times by* a factor of two or more. 
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The ADI procediire is a two-step, second-order accurate implicit integration 
technique. The first step advances the solution from the location (i) to an 
intermediate stage of computation (i+1/2) , much like a predictor step j, and 
the second step advances the solution to the new location (i+1) , much like a 
corrector step. The solution generated is then second order accurate about 
the point (i+1) . 

Two types of approaches are taken for obtaining the solutions of the governing 
equations (18) to (21) ; 

a. In the first method all four governing equations are solved simultaneously 
by the ADI method; 

b. In the second method three momentum equations [eqs. (18) to (20)] are 
solved simultaneously by the ADI method, and the Poisson-type pressure 
equation is solved by the Gauss-Sledel (explicit) method, 

A variable grid system is adopted for the finite difference formulation. It 
is shown in figure 3. The superscript i is- taken 'along axial direction x;-- the 
subscript j is taken along the n direction; and the subscript k is taken 
along :the ^ direction. The ratio of consecutive grid spacing along any axis 
is kept constant, i.e.. 




i 

- X 

AX 

i 

X - 

x^-l ■ 

An, 


- _ 

An_ 

- 

^i-1 


^k+l 

■ _ 

A?_ 

“ ^k - 

‘°k-l 


(31) 


In the variable grid system the derivatives of any dependent variable (j> 
are written in finite difference form as follows: 
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Figure 3. Discretized grid system and the ADI procedure. 
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9!i[ 

an2 


^n+ An_ 

2 — CAn^ + An_) 



l!i| 

3ti3 


j,k 


+. O(An^) 


( 33 ) 


A similar technique is used for the finite difference forms of and 

dQ 

9^ij) 

Linearization Scheme 

The ADI method is actually formulated for only linear equations. But the 
governing equations are nonlinear. Thus, a method must be developed to compute 
all the nonlinear coefficients (and the derivatives in the pressure equation) . 
This linearization can be accomplished by a quadratic extrapolation from the 
two previous x-stations i and (i-1). All of these nonlinear coefficients 
and derivatives are calculated at the (i+1/2) point by extrapolation using 



where Q refers to any of these quantities needed for extrapolation. 


19 



Boundary Conditions 


Boundary conditions o£ velocity components are very simple and are discussed 
in the previous section. Boundary conditions for pressure are taken as .or 
and are zero* In finite difference form to satisfy the above condition wall 

pressure px is denoted in terms of p 2 and p 3 by 


(1 + P2 - P3 

Pw = Pi = XCX + 2) 


(35) 


On the far field free boundaries the boundary conditions of all four variables 
are given during first iteration by equations (24) to (30) . After first iteration, 
the boimdary condition of all four variables is changed to 


3^(j) ^ 


For n direction 


and 


0 

9 ?^ 


For ? direction 


where (|) is one of the. variables u, v, w, p. Using finite difference formu- 
lation end value <!>„ is calculated in terms , and _ by 

^N-1 N-2 ^ 


(Pj^ = (1 + X)^ 


N- 


■ 1 ■ ^ ‘^N-2 


(36) 


Method (a) 

By using the extrapolation [eq. (33)] for the nonlinear coefficients and 
derivatives in equations (18) to (21) , the ADI difference approximation for 
these equations results in a linear set of four equations in four unknoTOS at 
each step of the marching procedure. A set of equations for the first half 
axial step and implicit in the n direction is given in Appendix B. A set 
of equations for the full axial step and implicit in the z direction is 
given in Appendix C. 
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For a more efficient approach to getting a solution, the four equations are 
coupled and solved simultaneously. The method is quite advantageous because, 
in addition to eliminating the choice of solution order, it makes more sense by 
allowing changes in any variable to be instantaneously sensed by all others. 

Let Wj^ = (u,v,w,p) be an unknown vector which contains four variables. For 
the sweep of the ADI procedure, where ^-derivatives are differenced implicitly 
along one line of the z direction, the resulting matrix equation derived from 
governing equations CIS) to (21) has a block tridiagonal structure as shown 
below. 


f^^k \-l \ l^^^k \+l 


(37) 


where 


= (u,v,w,p) 


(38) 


and Dj is a known source term vector along the n direction. 


Dj = (dll d2, d3, d4.)j 


(39) 


The coefficients (A), (B) , and (C) are block matrices of size 4x4, for 

example: 



an 

ax2 

ai3 

ai4 


azi 

a22 

^23 

^24 

[A],= 

&31 

^32 

as3 

^34 


r— 1 

1 

ait2 

a43 

^44 

The components 

a 

ran 

represent 


(40) 


equation at point k. The coefficients ^^^k given 

in Appendix D. 

The inversion of this block matrix is obtained by standard procedure . 
(ref. 13) . 
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Method (b] 


In this approach, instead o£ all four equations, three momentum equations 
are solved simultaneously with known values for pressure as calculated at station 
(i+1/2) by extrapolation of data from stations i and (i-1) . After the 
first iteration the new calculated pressure data of station (i+l/2) is used, 
obtained by interpolation of calculated pressure of station (i+1) and of 
known pressure of station i. 

During the first half axial step march equations (B-1) to CB-3^ from Appendix 
B are used, and during the next half axial step march the equations (C-1) to 
(C-3) from Appendix C are used. In this computation, the tridiagonal block 
matrices are of order 3x3 and the vector W = (u,v,w,). The coefficients of 
the matrices A, B, and C, and of source vector D are shoim in Appendix S. 

The pressure equation has no x derivative as the gradient of pressure 
along the axial direction is neglected. Thus, the pressure equation is an 
elliptic equation. Therefore, it is convenient to solve- the pressure equation 
by the Gauss-Siedel method. 


RESULTS AND DISCUSSION 

In the first phase of the computation, uniform space grid system (X = 1^ 
is adopted along the x, n, and 5 directions; 51 x 51 grid points are 
taken in the computation cross plane: h as well as ? are taken of 

the same magnitude and equal to 9. Thus, in the uniform space grid system, 
each grid space is of 0.18 x 0.18 size. The variables u, v, w, and p 
are calculated in the cross plane by providing the appropriate boundary condi- 
tions. In method (a), for better convergence, 5 to 10 iterations are adopted 
with the ADI method in each cross plane calculation. During the first iteration 
the coefficients of the nonlinear derivatives are extrapolated at the (i+1/2] 
location by the known values of i and locations, by equation (34) . 

Later -the nonlinear coefficients are calculated at (i+1/2) location by 
interpolation by the known values of i location and the latest computer values 
of (i+1) location. Near the leading edge 10 iterations are adopted, and in . 
the downstream part 5 iterations are considered. Some oscillations of the 
result are observed near the leading edge computation. In downstream marching, 
the oscillations are damped out. In the downstream section the result is observed 
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converged even in 5 iterations. In method C^) pressure equation (21) is solved 
by the Gauss-Siedel (explicit) method, thus for each iteration of ADI method 20 
iterations are given to the pressure equation. These are sufficient to obtain 
a converged solution. With both methods computation is carried out for a 
Reynolds number of 1000 and a Mach number 0.02. The results by both methods 
are almost identical. There is some deviation in pressure data-. The governing 
equations are completely independent of the value of the Reynolds ntimber. 

Figure 4 depicts isovels of the axial velocity u in the cross plane, no?. 
The results' found by Rubin and Grossman (ref, 4) are also plotted in the figure. 
Their result is obtained from similarity variables, and it is valid for any 
location along the axial direction. But in the present method the computation 
is carried out marching along the axial direction. The result plotted is for 
large values of axial distance. The isovels of the present result are very 
close to the isovels of Rubin and Grossman. The present result is in good agree- 
ment with their result. The result indicates that the boundary layer thickness 
6 on the planes n = 0 and ? = 0 toward the corner layer and both the perpen- 
dicular boundary layers merge to a single boxmdary layer by a smooth curvature. 
More details about the characteristics of the axial velocity are discussed in 
references 2 to 5. Thus, the discussion about the axial velocity is not given 
here. 

Axial velocity u along the symmetry line is plotted in figure 5. The, plot 
reveals that in the corner layer, due to the growth of the crossflow, there is 
a decrement of the axial velocity. • 

The skin friction parameter at n = 0 is plotted along the ? direction 

in figure 6. The skin friction coefficient is defined by 

_ / 2v \ 

^f “ yujc J 3n'n=0 

as shown in figure 6. The skin friction parameter is zero at the corner and 
increases monotonically to its asymptotic two-dimensional value of 0.4696. 

This result is compared with the result of Rubin and Grossman, and the present 
result ~is in good agreement with their result. 

The lateral velocity w along the symmetry line is plotted in figure 7. 

The velocity w increases rapidly in the corner region to a maximum value at 
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Figure 4. Isovels o£ the axial velocity. 
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the edge of the boxmdary layer. In the potential flow region this velocity decays 
algebraically toward the value of zero.- It is observed that cross velocities 
near the leading edge are very small. The crossflow velocities^ however, increase 
with increased axial distance, i.e., a secondary flow develops. This crossflow 
effect is of a significant amount (although constant) in the downstream portion 
of the comer. 

The secondary crossflow vectors are plotted in figure 8. The traces inside 
the corner layer are drawn towards' the comer line. These suggest the presence 
of the transverse flow along the comer walls which is towards the corner line; 
i.e., outflow from the corner. This plot also predicts the presence of vortices 
on both sides of the symmetry line, but the vortices are zero along the symmetry 
line. 

The pressure contours are plotted in figure 9. Within the corner layer 
pressure is almost constant. Near the boundary layer edge the pressure decreases 
rapidly and reaches the minimum value. Then again, pressure rises slightly in 
the potential flow region. The results obtained by method (a) and method (b) 
are compared in this figirre. The results are not in complete agreement at all 
grid points . 

In figure 10 the effect of the axial distance on axial velocity is shoivn. 

Tlie axial velocity along the symmetry line is plotted at different values of x 
location. The changes in the axial velocities are noted as the axial distance 
X is increased. This plot indicates that similarity solution is valid only 
for large values of axial distance. 

CONCLUDING REMARKS 

The flow along a comer formed by the intersection of two perpendicular flat 
plates leads to a three-dimensional boundary layer problem. 

It is concluded that the similarity solution is indeed valid, but only in 
the dovmstream section. This solution, however, is not valid near the leading 
edge. 
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Figure 8. Secondary flow vectors in the comer layer. 
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APPENDIX A 


POTENTIAL FLOW SOLUTION 


After Ghia and Davis (ref, 5) the three velocity components in dimensionTess 
form are given by: 


u = 1 


2v^ 


y{x^ - y^} 


- 1/2 


I X I 


1/2 


2v^ 


:(x^ + z^) 


- 1/2 


I X +Vx^ + I 


1/2 


(A-1) 


V = v(x,y) = 


2^ 


-t-Vx^ + y^ 
x^ + y^ 


1/2 


(A-2) 


W = Vf(x,z3 = 




X + + Z'^ 

x^ + z^ 


1/2 


CA-3) 


From (x,y, z) coordinate system u, V, and ^ are transformed to the new 
coordinate (xjH,?) system, and the transformed velocity components, u, v, and 
w, respectively are obtained as: 
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u = 1 


3 

Re ^ 



(A-4) 


g 

Re 




1/2 


V = v(x,n) 


_g_ 

/2 


X 



1/2 


X + 


2n^ 

Re 


(A-5) 


w = w(x,? ) = — 


e. 


X + 


I 


X + 


Re 


2c^x 

Re 


1/2 


(A-6) 


The results obtained above are the same as those determined by Rubin 
(ref. 2) by using Green's function. 
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APPENDIX B 


HALF AXIAL STEP - IMPLICIT IN n DIRECTION 


o... 9u* ^ ^ 9u* , _ , 9u^ 

2xu -Trr- - (un - v) — - (u? - w) = + 

-5^2 3j,2 


3x 




CB-1) 


2xu 


9v* 

9X 


, ^ 9v* , ^ , 9v 

Cu„ - V) -5^ - C15 - W) -5^ 


, 3p» S^v* 

• 2x “T — + + 

9n2 3?^ 


UV*^ 


CB-2) 


^ 9w* , ^ ^ ^ 9w ^ 

-9T - -9TT ■ TT ■ 


= -2x ^ + 
^ 9C 


9^W^ + 9^w^ 

,2 


9n' 


dZ“ 


(B-3) 


9^p* . 9^p 1 

» " ■ ' J5JJ 


3n" 


9?^ 


X 


„ 9u* 9u* 

2x -;rr- - n 


9x 


9p 


9u 

^ Itr 


9v 9w 
9p 9? 


, 9v 9w 3u 9v* 9v* , 9v . 

9n 9? ~ 9n 9x ■ ^ 9 ti ' ^ 9? ~ ^ 


CB-4) 


9v 9w 9u 9w* 9w* ^ 9w , 

9? 9t) " 9? 9x ■ ^ 9n ■ ^ 9^ ' 


where the terms shown without superscripts are considered at the station (i+1/2) 
The variables shown by superscript * are considered unknown at station (i+1/2) 
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APPENDIX C 


FULL AXIAL STEP - IMPLICIT IN t, DIRECTION 


2x u 


~ 1+1 ~ ^ 1+1 
3u f . > 3u* , ^ . 3u 

- (un - v) - (u? - w) 


3 X 


3n 


3? 


3^u* ^ 3^u 


3n" 


2,-i+l 


3?" 


(C-ID 


2x u 


- Cun - v) - (uc - w) 


3x 


3n 


3? 


- UV* 


= -2x 


9e1 + 

3n 



3?^ 


(C-2) 


2X u 


3w 


i+1 


3x 


- Cun - w) 


3w* 

3n 


Cus - w) 


3w 


i+1 


3? 


- uw*' 


= -2x 


3£ 


i+1 


3? 


2-tct'* 3^w 

+ — 


3^w 


i-rl 


9ti^ 




(C-3) 


3^p* ^ 

3n^ 35^ 



n 


3u* 

3n 




3u^\ /3v 
3C / \3n 



3w 3u f 3 v^~^^ 3v* 

3p 3? ■ 3n 3x ■ 3n 



3v 3w 3u /- 3w^^^ 

3? 3n ' 3C 3x 


3w* 

"1^- 




3w 


i+1 


3? 



CC-4) 
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where the terras shotTO without superscripts are considered at the station 
The variables shown by the superscript are considered known at station 
(i+1/2) . The variables shown by the superscript Ci+i) 3^^® only the unknown 
variables at station (i+l) . 
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APPENDIX D 


MATRIX COEFFICIENTS OF THE LINEARIZED DIFFERENCE EQUATIONS 

METHOD Ca) 


For the first half of the ADI procedure, from a location i to the inter- 
mediate * location, the components of matrices (A), (B) , (C) , and the 

vector of the set of equations of Appendix B are as follows: . 


an 


(up - 
■ CX + 


2X . 

(X + i)An^ 
+ 


CD-I) 


ai2 = 0 


CD-2)- 


ai3 = 0 


CD-3) 


ai4 = 0 


CD-4) 


a2i - 0 


a22 = aix 


a23 = 0 


2X^x 

■ cx + l)An^ 


CD-5) 


CD-6) 


CD-7) 


CD-8) 
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asi = 0 


CD-9) 


a32 = 0 


CD-10) 


^33 = ail 


CD-11) 


as4 = 0 


CD-12) 


a4i 




3v ^ 3w 




3n 3? I CA + l)An^ 


. - 1 3 u I nx^ \ 

X 3ri \ CX + 1)Ati^ j 


^4 3 


1 3u i nX^ \ 
“x 3? (CX + 1)Ati I 


a44 — 


2X 


CX + i)An2 


CD- 13) 


(D-14) 


CD-15) 


CD-16) 


bn 


4xu 

Ax 


(un - v) CX - 1) 
Ad. 



CD- 17) 
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bl2 = 0 CD-18) 

bi3 = 0 (D-19) 

bi4 = 0 (D-20) 

b2i = 0 CD-21) 

b22 = bii - CD-22) 

b23 = 0 CD-23) 

hik = CD-24) 

+ 

b3i = 0 CD-25) 

b32 = 0 (D-26) 
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bs3 = i>Z2 


CD-27) 


bs4 - 


(D-28) 



lix. + Sw W 4x 
|3n I \ Ax 


n(A - 1) \ 

( 


CD-29) 


^42 


= illl /4x 
X 3n 1 Ax 


nCA - 1) 

An^ 


- 1 


CD-30) 


b43 


X 3? ( 


■4x 

Ax 


nCX - 1) 

Ad, 


CD-31) 


b44 = 



CD-32) 


Cii = - 


Cup - v) 

(A + l)An^ 


. 2 

CA + i)An^ 


CD-33) 


Ci2 = 0 


CD-34) 


C13 = 


CD- 35) 
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Cl4 = 0 


CD-36) 


C21 = 0 


CD-37) 


C 22 = Cii 


CD-38) 


C23 = 0 


CD- 39) 


C24 = 


2x 


CA + l)An^ 


CD-40) 


C31 = 0 


CD-41) 


C 32 = 0 


CD-42) 


C33 == cii 


CD-43) 


C34 = 0 


CD-44) 
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CD-45) 


C 41 


C 42 


C43 


C44 


di = 


d2 = 


da = 


d4 = 


1. n 

X 1 3n 3<; i (X + i)Ati^ 


^ 3ju ri 

X 9n c^ 


(D-46) 


^ n 

X H (X + l)An_j_ 


(D-47) 


2 

(X + i)AnJ 


(D-48) 


4xuu^ ^ 
Ax 




4xuv^ 

Ax 


+ (u? 


4Xuw^ 

Ax 


+ (u? 


w) 


w) 


w) 


3u^ 

i 


9C 

9?2 

3v^ 

3^v^ 

^ 'T' 

9? 

9C^ 

3w^ 


9C 

9C^ 


- 2x 


3£i 

9? 


(D-49) 


(D-50) 


(D-51) 


9^p^ ^ £ 3u^ 3V 3w ^ 3v 3w 

gj .2 X 9? \3n 3? j ■ 9n 3c 3?: 3n 

(D-52) 

C 9u 3v^ ^ 3u 3w^ 

X 3n 3? * X ^r, 3? 
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where the unsubscripted field variables are evaluated at the point (i+l/2,j,k) 
by the extrapolation equation (34] during first iteration and by the interpola- 
tion after one iteration using the calculated riew value of previous iteration. 

For the second half of the ADI procedure, from the location * to the 
(i+1) location, the same standard procedure is adopted. Here the derivatives 
of T) axe denoted by *, and they are assumed known, and the derivatives of 
Z are assumed at station (i+1) and are asstmied unknown. Thus, for the second 
half of the ADI procedure the elements of the matrices (A) , (B) , (C) , and 

the vector are. not given. 
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APPENDIX E 


MATRIX COEFFICIENTS OF THE LINEARIZED DIFFERENCE EQUATIONS 

METHOD (b) 


As discussed in the Numerical Analysis Section, in this approach the 
equations (18) to (20) are solved by the ADI method assuming the known pressure 
values. Thus, the three equations make the block matrices of size 3x3. The 
elements of the block matrices (A) , (B) , and (C) are identical to those 
of Appendix D. The constant vector 6^ is slightly different, and the elements 
of vector D are given below. 


j 4xuu' 

dj = 


dz 


AX 


do = 


4xuv 



3^u^ 


3?2 


32v^ 

3? 

3C^ 

3w^ , 

*r 



3n 


2x 




Ax 3C • 3^2 3? 


(E-1) 


(E-2) 


(E-3) 


A similar approach is adopted during the second half of the ADI procedure. During 

both, the procedures and are computed at (i+l/2,j,k). 

3ti 0? 

The pressure is calculated from equation (21) after obtaining the solution 
of u, V, and w from equations (18) to (20) by the explicit method. The 
Gauss-Siedel method is used for solution of the pressure equation. The complete 
expression for pressure calculation is as follows: 
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9v 3w 9u 

3c 3n ~ 9^ 
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